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abstract 

We present the 1667 MHz OH megamaser luminosity function derived from a single flux-limited sur- 
vey. The Arecibo Observatory OH megamaser (OHM) survey has doubled the number of known OH 
megamasers, and we list the complete catalog of OHMs detected by the survey here, including three rede- 
tections of known OHMs. OHMs are produced in major galaxy mergers which are (ultra)luminous in the 
far-infrared. The OH luminosity function follows a power law in integrated line luminosity, <f> oc L^jf^ 



Mpc~ dex~ , and is well-sampled for 10 < Lqh < 10 L©. The OH luminosity function is incor- 
porated into predictions of the detectability and areal density of OHMs in high redshift OH surveys for 
a variety of current and planned telescopes and merging evolution scenarios parameterized by (1 -f z)™ 
in the merger rate ranging from m = (no evolution) to m = 8 (extreme evolution). Up to dozens of 
OHMs may be detected per square degree per 50 MHz by a survey reaching an rms noise of 100 /iJy 
per 0.1 MHz channel. An adequately sensitive "OH Deep Field" would significantly constrain the evo- 
lution exponent m even if no detections are made. In addition to serving as luminous tracers of massive 
mergers, OHMs may also trace highly obscured nuclear starburst activity and the formation of binary 
supermassive black holes. 

Subject headings: masers — galaxies: interactions — galaxies: evolution — galaxies: luminosity 
function — galaxies: starburst — radio lines: galaxies 



1. INTRODUCTION 

OH megamasers are luminous masing lines at 1667 and 
1665 MHz which are at least a million times more lumi- 
nous than typical OH masers associated with compact HII 
regions. OH megamasers (OHMs) are produced in (ul- 
tra)luminous infrared galaxies ([UjLIRGs), major galaxy 
mergers undergoing extreme bursts of circumnuclear star 
formation. OH megamasers are especially promising as 
tracers of dust obscured star formation and merging galax- 
ies because they can potentially be observed at high red- 
shifts with modern radio telescopes in reasonable integra- 
tion times, they favor regions of high dust opacity where 
ultraviolet, optical, and even near-IR emission can be ex- 
tremely difficult to detect, and their detection automat- 
ically provides an accurate redshift measurement. The 
Arecibo Observatory^ OH megamaser survey is a flux- 
limited survey designed to quantify the relationships be- 
tween merging galaxies and the OH megamasers (OHMs) 
they produce with the goal of using OHMs as luminous 
tracers of mergers at high rcdshifts (Darling & Giovanelli 
2000, 2001, 2002; hereafter Papers I, II, III). Central to 
the application of OHMs as tracers of merging galaxies at 
various redshifts is a measurement of the low redshift OH 
luminosity function. 

The Arecibo OH megamaser survey is the first sur- 
vey with adequate statistics to construct an OH lumi- 
nosity function (LF) from a flux-limited sample. Baan 
(1991) computed a OH luminosity function from the 48 
OHMs which were known at the time under the assump- 
tion that they constituted a flux-limited sample and that 
the fraction of LIRGs showing detectable OHMs is con- 
stant. Baan's OH LF has a Schecter-type profile which 
is slowly falling from Lqh — 10°-10^ Lq, has a knee at 

^ The Arecibo Observatory is part of the National Astronomy and 
Ionosphere Center, which is operated by Cornell University under a 
cooperative agreement with the National Science Foundation. 



roughly 10^-^ Lq, and a steep falloff out to 10^ Lq. We 
are now in a position to recompute the OH LF from a 
single complete survey with well-defined selection criteria 
extracted from the PSGz, a flux-limited catalog which also 
has well-defined selection criteria (Saunders et al. 2000). 
An OH LF will be a useful guide for deep surveys for OHMs 
which can be related to the merging history of galaxies, 
the dust-obscured star formation history of the Universe, 
and the production of some portions of the low frequency 
gravitational wave background. 

This paper presents an overview of the Arecibo OH 
megamaser survey, focusing on the issues pertinent to con- 
structing an OH luminosity function (OH LF) from the 
survey results. Selection methods and the complete cata- 
log of detected OHMs are presented in §2. Section 3 dis- 
cusses the methods used to compute the OH luminosity 
function, presents the OH LF, and compares it to previ- 
ous resuhs for OHMs and ULIRGs. The OH LF is then 
applied to the problem of detecting OHMs at high redshift 
in §4 and some discussion is made of the utility of OHMs 
as tracers of galaxy evolution, including merging, dust- 
obscured nuclear starbursts, and the formation of binary 
supermassive black holes. 

Note that Papers I-III assume a cosmology with Ho = 
75 km s^^ Mpc^^, qo — 0, and Ha = for ease of com- 
parison with previous (U)LIRG surveys such as Kim & 
Sanders (1998). This analysis, however, assumes a more 
likely cosmology which is fiat but accelerating: S7m = 0.3 
and — 0.7. All of the OHM data presented here have 
been converted to this cosmology. 

2. THE ARECIBO OH MEGAMASER SURVEY 

Papers I-III present a survey for OHMs conducted at 
the Arecibo Observatory covering one quarter of the sky 
to a depth of roughly one gigaparsec. The survey doubles 
the number of known OH megamasers and quantifies the 
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relationship between luminous infrared galaxies and OH 
megamasers. The survey builds the foundation required to 
employ OH megamasers as tracers of major galaxy merg- 
ers, dust-obscured star formation, and the formation of 
binary supermassive black holes spanning the epoch of 
galaxy evolution to the present. Here we present details of 
the survey relevant to constructing an OH LF: the candi- 
date selection criteria and the complete catalog of detected 
OHMs with properties converted to a flat Cl\ = 0.7 cos- 
mology. 

2.1. Candidate Selection 

For the Arecibo OH megamascr survey, candidates 
were selected from the Point Source Catalog redshift sur- 
vey (PSCz; Saunders et al. 2000), supplemented by the 
NASA/IPAC Extragalactic Database.^ The PSCz catalog 
is a flux-limited (IRAS /eo/^m > 0.6 .Jy) redshift survey of 
15,000 IRAS galaxies over 84% of the sky (sec Saimders 
et al. 2000). We select IRAS sources which are in the 
Arecibo sky (0° < ^ < 37°), were detected at 60 /zm, and 
have 0.1 < z < 0.45. The lower redshift boimd is set to 
avoid local radio frequency interference (RFl), while the 
upper bound is set by the bandpass of the wide L-band 
receiver at Arecibo, although an effective upper bound 
is imposed around z = 0.23 by the radio frequency in- 
terference (RFI) environment, as discussed in §3.3. No 
constraints are placed on FIR colors or luminosity. The 
redshift requirement limits the number of candidates in 
the Arecibo sky to 311. The condition that candidates 
have z > 0.1 automatically selects (U)LIRGs if they are 
included in the PSCz. The strong influence of LpiR on 
OHM fraction in LIRGs (see Paper HI) is the primary 
reason for our high detection rate compared to previous 
surveys (e.g. Staveley-Smith et al. 1992; Baan, Haschick, 
& Henkel 1992). 

2.2. OH Megamaser Detections 

Tables 1 and 2 list respectively the optical/FIR and ra- 
dio properties of the 50 new OHM detections and 3 re- 
detections. Note that OHM IRAS F11180+1623 is not 
in the PSCz sample, but was observed along with other 
OHM candidates not found in the PSCz sample to flll in 
telescope time when local sidereal time coverage of the of- 
flcial sample was sparse. This detection is not included 
in any of the survey statistics or interpretation, including 
the OH LF. Spectra of the 53 OHMs appear in Papers 
I, II, and HI. Table 1 lists the optical redshifts and FIR 
properties of the non-detections in the following format: 
Column (1): IRASY&vat Source Catalog (FSC) name. 
Columns (2) and (3): Source coordinates (epoch B1950.0) 
from the FSC, or the Point Source Catalog (PSC) if un- 
available in the FSC. 

Columns (4), (5) and (6): Heliocentric optical redshift, 
reference, and corresponding velocity. Uncertainties in ve- 
locities are listed whenever they are available. 
Column (7): Cosmic microwave background rest-frame ve- 
locity. This is computed from the heliocentric velocity us- 
ing the solar motion with respect to the CMB measured 

2 The NASA/IPAC Extragalactic Database (NED) is operated by 
the Jet Propulsion Laboratory, California Institute of Technology, 
under contract with the National Aeronautics and Space Adminis- 
tration. 



by Lineweaver et al. (1996): czq = 368.7 ± 2.5 km s ^ 
towards (Lb) = (264?31 ± 0?16,48?05 ± 0?09). 
Column (8): Luminosity distance computed from zcmb 
via 

„ rZCMB _i 

Dl = {1 + zcMB)-rr / [(1 + z'f^M + Oa] ' dz', 

tlo Jo 

(1) 

assuming flM = 0.3 and 0,a = 0.7. 

Columns (9) and (10): IRAS 60 and 100 /im flux densi- 
ties in Jy. FSC flux densities are listed whenever they are 
available. Otherwise, PSC flux densities are used. Un- 
certainties refer to the last digits of each measure, and 
upper limits on 100 /im flux densities are indicated by a 
"less-than" symbol. 

Column (11): The logarithm of the far-infrared lumi- 
nosity in units of Hj^Lq. Lpm is computed following 
the prescription of Fullmer & Lonsdale (1989): Lfir= 
3.96 X 10'"^ £'1(2.58/60 + /loo), where /eo and /loo are the 
60 and 100 /Ltm flux densities expressed in Jy, Dl is in 
h^^ Mpc, and Lpm is in units of H^^Lq. If /loo is only 
available as an upper limit, the permitted range of LpiR is 
listed. The lower bound on LpiR is computed for /loo = 
Jy, and the upper bound is computed with /loo set equal to 
its upper limit. The uncertainties in and in the IRAS 
flux densities typically produce an uncertainty in log LpiR 
of 0.03. 

Table 2 lists the OH emission properties and 1.4 GHz 
flux density of the OH detections in the following format: 
Column (1): IRAS FSC name. 

Column (2): Measured heliocentric velocity of the 1667.359 

MHz line, defined by the centcT of the FWHM of the line. 
The uncertainty in the velocity of the line center is esti- 
mated assuming an uncertainty of ±1 channel (±49 kHz) 
on each side of the line. Although one can generally de- 
termine emission line centers with much higher precision 
when the shapes of lines are known, OHM line profiles 
are asymmetric, multi-component, and non-gaussian (see 
Papers I-III). They defy simple shape descriptions, so we 
use this conservative and basic prescription to quantify 
the uncertainty in the line centers. 

Column (3): Peak flux density of the 1667 MHz OH fine 
in mjy. 

Column (4): Equivalent width-like measure in MHz. 
M^i667 is the ratio of the integrated 1667 MHz line flux 
to its peak flux. Ranges are listed for VKiee? in cases 
where the identification of the 1665 MHz line is unclear, 
but in many cases the entire emission structure is included 
in PFi667 as indicated in the discussion of each source in 
Papers I~III. 

Column (5): Observed FWHM of the 1667 MHz OH line 
in MHz. 

Column (6): Rest frame FWHM of the 1667 MHz OH line 
in km s~^. The rest frame width was calculated from the 
observed width as Av^est = c(l + z){Avobs/t^o)- 
Column (7): Hyperfine ratio, defined by Rh — -Fi667/-Pi665) 
where iv is the integrated flux density across the emission 
line centered on v. Rh = 1-8 in thermodynamic equi- 
librium. In many cases, the 1665 MHz OH line is not 
apparent, or is blended into the 1667 MHz OH line, and 
a good measure of Rh becomes difficult without a model 
for the line profile. It is also not clear that the two lines 
shotild have similar profiles, particularly if the lines are 
aggregates of many emission regions in different satura- 
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Table 1 

Arecibo Survey OH Megamasers: Optical Redshifts and FIR Properties 
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''A known OHM included in the survey sample. 

^ IRAS 11180+1623 is not in the PSCz catalog (excluded by the PSCz mask; Saunders et al. 2000). It is in the 1 Jy survey (Kim & Sanders 
1998). 

References. — Redshifts were obtained from: (1) Saunders et al. 2000; (2) Strauss et al. 1992; (3) Lawrence et al. 1999; (4) Lu & Freudling 
1995; (5) Kim & Sanders 1998; (6) Fisher et al. 1995. 
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Table 2 



Arecibo Survey OH Megamasers: OH Line and 1.4 GHz Continuum Properties 



IRAS Name 




/l667 


Wl667 






Rh 


log Lj?/ij 


log Lojr 


fl.4GHz 


FSC 


km/s 


mJy 


MHz 


MHz 


km/s 








mJy 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 



01562+2528 


49814 


15) 


6.95 


1.29 


1.04 


218 


5.9 


11.90 


3.25 


6.3(0.5) 


02524+2046 


54162 


15) 


39.82 


0.50 


0.36 


76 


3.2 


11.81-12.28 


3.74 


2.9(0.6) 


03521+0028 


45512 


14) 


2.77 


0.61 


0.29 


59 


5.8 


12.28 


2.44 


6.7(0.6) 


03566+1647 


39865 


14) 


1.96 


0.98 


0.23 


48 


> 9.6 


11.40-11.75 


2.32 


3.5(0.5) 


04121+0223 


36590 


14) 


2.52 


0.76 


1.04 


209 


2.9 


11.40-11.69 


2.32 


3.1(0.5) 


06487+2208 


42972 


14) 


6.98 


0.66 


0.37 


76 


8.0 


12.10 


2.82 


10.8(0.6) 


07163+0817 


33150 


14) 


4.00 


0.69 


0.12 


24 


~ 5.5 


11.53 


2.37 


3.5(0.5) 


07572+0533 


56845 


15) 


2.26 


1.03 


0.73 


156 


> 10.4 


12.04 


2.74 


< 5.0 


08201+2801 


50325 


15) 


14.67 


0.97-1.19 


0.98 


205 


> 8.2 


12.00 


3.45 


16.7(0.7) 


08279+0956 


62422 


15) 


4.79 


1.02 


0.95 


207 


5.9 


11.75-12.01 


3.23 


4.4(0.8) 


08449+2332 


45424 


14) 


2.49 


1.09 


0.47 


97 


> 11.0 


11.79 


2.59 


6.1(0.5) 


08474+1813 


43750 


14) 


2.20 


1.29-1.70 


1.98 


409 


> 3.0 


11.91 


2.70 


4.2(0.5) 


09039+0503 


37720 


14) 


5.17 


1.23 


1.05 


212 


8.5 


11.86 


2.83 


6.6(0.5) 


09531+1430 


64434 


15) 


3.98 


1.03 


1.17 


256 


~ 3.4 


12.08 


3.42 


3.0(0.5) 


09539+0857 


38455 


14) 


14.32 


1.47 


1.56 


317 


2.5 


11.79 


3.48 


9.5(1.2) 


10035+2740 


50065 


14) 


2.29 


0.74 


0.31 


65 


15.4 


12.01 


2.50 


6.3(0.5) 


10339+1548 


58983 


15) 


6.26 


0.28 


0.19 


40 


^ 14.5 


12.09 


2.65 


5.1(0.5) 


10378+1108'' 


40811 


14) 


19.70 


1.50 


0.87 


177 




12.05 


3.54 


8.9(0.6) 


11028+3130 


59619 


15) 


4.27 


0.72 


0.41 


89 


5.5 


12.13 


2.97 


< 5.0 


11180+1623" 


49783 


14) 


1.82 


0.42 


0.61 


127 


> 5.1 


12.02 


2.34 


4.2(0.5) 


11524+1058 


53404 


15) 


3.17 


1.21 


1.32 


279 


~ 4.9 


11.93 


2.98 


< 5.0 


12005+0009 


36472 


14) 


3.51 


0.71 


0.41 


82 


~ 2.0 


11.52 


2.62 


5.4(0.6) 


12018+1941'' 


50317 


15) 


2.63 


0.81 


0.86 


181 


> 5.6 


12.16 


2.59 


6.5(0.5) 


12032+1707 


64920 


15) 


16.27 


2.69 


3.90 


853 




12.31 


4.15 


28.7(1.0) 


12162+1047 


43757 


141 


2.07 


0.64 


0.51 


105 


> 11.1 


11.50-11.68 


2.25 


6.8(1.6) 


12549+2403 


39603 


14) 


1.79 


0.89 


0.50 


102 


^ 2.6 


11.60 


2.37 


3.7(0.5) 


13218+0552 


61268 


15) 


4.01 


2.49 


1.45 


314 




12.13 


3.45 


5.3(0.5) 


14043+0624 


33912 


14) 


2.75 


0.33 


0.27 


54 


1.4 


11.51 


2.10 


15.6(1.0) 


14059+2000 


37246 


14) 


15.20 


1.10 


0.80 


161 


5.3 


11.68 


3.34 


7.5(0.5) 


14070+0525'* 


79929 


16) 


8.37 


3.21 


2.62 


596 




12.55 


4.13 


4.0(0.6) 


14553+1245 


37462 


14) 


2.93 


0.39 


0.38 


77 


> 14.5 


11.62 


2.27 


3.8(0.5) 


14586+1432 


44380 


14) 


7.11 


< 2.67 


1.79 


369 




11.64 


3.41 


11.1(0.6) 


15224+1033 


40290 


14) 


12.27 


0.73-0.80 


0.15 


31 


> 9.5 


11.57 


3.04 


3.6(0.5) 


15587+1609 


40938 


14) 


13.91 


0.99 


0.86 


176 


6.9 


11.60 


3.26 


< 5.0 


16100+2527 


40040 


14) 


2.37 


0.60 


0.23 


46 


3.2 


11.39-11.63 


2.29 


< 5.0 


16255+2801 


40076 


14) 


7.02 


0.45 


0.39 


79 


> 13.7 


11.69 


2.57 


< 5.0 


16300+1558 


72528 


15) 


3.12 


0.56 


0.59 


131 




12.47 


2.85 


7.9(0.5) 


17161+2006 


32762 


14) 


4.84 


0.62 


0.38 


76 


~ 6.2 


11.16-11.43 


2.39 


7.3(0.6) 


17539+2935 


32522 


14) 


0.76 


0.72 


0.81 


161 


> 2.9 


11.58 


1.76 


4.0(0.6) 


18368+3549 


34832 


14) 


4.58 


1.79 


2.10 


421 


~ 9.5 


11.98 


2.85 


21.0(0.8) 


18588+3517 


31686 


14) 


7.37 


0.56 


0.32 


64 


5.1 


11.66 


2.52 


5.9(0.5) 


20248+1734 


36538 


14) 


2.61 


1.36 


0.88 


177 


~ 6.8 


11.68 


2.53 


< 5.0 


20286+1846 


40471 


14) 


15.58 


1.51 


1.10 


224 


> 4.4 


11.81 


3.41 


< 5.0 


20450+2140 


38398 


14) 


2.27 


0.67 


0.71 


144 


> 6.2 


11.67 


2.24 


5.0(0.5) 


21077+3358 


52987 


15) 


5.04 


1.86 


1.15 


243 


> 7.4 


11.75-11.98 


3.26 


9.4(1.0) 


21272+2514 


45032 


14) 


16.33 


1.87 


1.27 


263 


13.7 


11.69-11.89 


3.66 


4.4(0.5) 


22055+3024 


37965 


14) 


6.35 


0.77 


0.46 


92 


6.2 


11.93 


2.73 


6.4(0.5) 


22116+0437 


58180 


15) 


1.76 


1.16 


0.56 


121 


~ 5.2 


11.85-12.01 


2.77 


8.4(0.6) 


23019+3405 


32294 


14) 


3.58 


0.52 


0.28 


57 


> 15.6 


11.69 


2.12 


7.7(0.5) 


23028+0725 


44529 


14) 


8.69 


1.09 


1.06 


219 


1.9 


11.61-11.81 


3.29 


19.5(1.1) 


23129+2548 


53394 


15) 


4.59 


2.0 


1.78 


376 




12.21 


3.27 


4.7(0.5) 


23199+0123 


40680 


14) 


1.80 


0.82 


0.68 


139 


2.3 


11.58 


2.38 


3.0(0.5) 


23234+0946 


38240 


14) 


3.32 


1.23 


1.32 


266 


2.4 


11.88 


2.75 


11.6(1.0) 



*The frequency width Aviqqj is the observed FWHM. 

The velocity width Avi667 is the rest frame FWHM. The rest frame and observed widths are related by Avrest — 

c(l + z){Aiy,ts'/i^o). 

^1.4 GHz continuum fluxes are courtesy of the NRAO VLA Sky Survey (Condon et al. 1998). 

■^Redetection of a known OHM included in the survey sample. 

^IRAS 11180+1623 is not in the PSCz catalog (excluded by the PSCz mask; Saunders et al. 2000). It is in the 1 Jy survey 
(Kim & Sanders 1998). 
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tion states. Some spectra allow a lower limit to be placed 
on Rh, indicated by a "greater than" symbol. Blended or 
noisy lines have uncertain values of Rh, and are indicated 
by a tilde, but in some cases, separation of the two OH 
lines is impossible and no value is listed for Rh- 
Column (8): Logarithm of the FIR luminosity, as in Table 
1. 

Column (9): Logarithm of the measured isotropic OH line 
luminosity, which includes the integrated flux density of 
both the 1667.359 and the 1665.4018 MHz lines. 
Column (10): 1.4 GHz continuum fluxes from the NRAO 
VLA Sky Survey (Condon et al. 1998). If no continuum 
source lies within 30" of the IRAS coordinates, an upper 
limit of 5.0 mJy is listed. 

3. THE OH LUMINOSITY FUNCTION 

3.1. Computing a Spectral Line Luminosity Function 

The luminosity function i>(i) is the number density of 
objects with luminosity L per (logarithmic) interval in L. 
An unbiased direct measurement of $(i) would require 
that all objects with a given luminosity be detected within 
the survey volume, which is generally not possible in a 
flux-limited survey. Instead, each object in a survey has 
an effective volume in which it could have been detected 
by the survey, and the sum of detections weighted by their 
available volumes Va determines the luminosity function. 
The most general unbiased maximum likelihood computa- 
tion of a luminosity function is computed from Va following 
the prescription (Page & Carrera 2000): 

where the sum is over the detected objects in the luminos- 
ity bin A log L centered on L and redshift bin Az centered 
on z. The uncertainty in the luminosity function is 



1 



Computation of the volume available to each detection de- 
pends on the areal coverage of the survey f2 and the maxi- 
mum detectable distance of each object detected {DL,max 
or z.fnax)- The relationship between physical comoving vol- 
ume and luminosity distance is given by 

3 



V = 



1 + z 



(4) 



only when flk = 0, where fife = 1 — Hm — (Weinberg 
1972). From the definition of luminosity distance. 



Dr 



\4ttS 



(5) 



for fife = 0, where L is the integrated line luminosity and 
S is the integrated line flux. Detection of spectral fines 
depends on the peak line flux density rather than the in- 
tegrated line flux. The luminosity distance of Equation 5 
must be modified by a factor 1 + z in order to change to 
rest-frame luminosity density and flux density: 



Dr 



{i + zY- 



(6) 



For a survey with sensitivity limit na^, where is the 
RMS noise flux density at frequency and assuming that 



this noise is independent of frequency, L^^ is an invariant 
quantity which relates the observed luminosity distance 
and peak flux density to the maximum detectable dis- 
tance and the sensitivity limit ncr^. We thus obtain an 
expression for the maximum detectable distance, DL^maxi 
which depends on itself through z^ax 



D 



L,max 



= Dt 



l + z 



(7) 



and must be solved such that D^^max can be related to ob- 
served quantities only. Equation 1 has no simple analytic 
solution for Cl\ ^ and must be solved numerically. In- 
serting this solution into Equation 7, we obtain a relation 



m Z^ax' 



/ [(! + • 
Jo 



z' f^M + ^A\ ' dz' 



(8) 



= — === {l+Zmax) ^ 

cVl + z \na^ J 
which has no analytical solution. The numerical solution 
for Zjnax (and the equivalent fi^mox) can be inserted into 
an expression for Vmax obtained from Equation 4. Finally, 
we obtain the volume available to a given emission line 
source in the survey in the case where there is a minimum 
survey redshift Zmin'- 

Va — Vmax Vmin 

3 / \ 3" 



3" 



l + z„ 



1 + Zr, 



(9) 



Note that if Zmax is greater than the upper bound in red- 
shift of the survey, then we set Zmax equal to this upper 
bound. 

3.2. Computing a Continuum Luminosity Function 

The procedure for computing a continuum luminosity 
function closely follows the steps for the emission line com- 
putation. The detection threshold is determined in the 
same manner, from a flux density sensitivity rather than 
an integrated flux. Unlike spectral line surveys, however, 
continuum surveys do not usually tune the bandpass to 
the redshift of each observed object. This introduces an 
error in computing the invariant luminosity from the mea- 
sured flux because a different portion of the rest-frame 
spectral energy distribution is sampled for each object 
due to the redshift. A k-correction must be applied to 
the measured flux density to obtain the correct luminosity 
and hence to compute the maximum detectable distance 
for each source. The net k-correction between z and Zmax 
will modify Equation 7 slightly 



Dl 



1 



l + z 



(10) 



but will add another wrinkle to the derivation of V^ax 
in §3.1 because the k-correction will itself depend on Dl. 
We will approximate the k-correction to be constant at 
first, determine D^ max, recompute the k-correction, and 
then correct DL,max if necessary. The k-correction is a 
weak function of Dl — DL.max, and this single iteration 
approach to computing DL,max is adequate. 
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3.3. The OH Megamaser Luminosity Function 

The Arecibo OH megamaser survey is a flux-limited sur- 
vey for OH emission lines, but has three additional con- 
straints: (1) objects observed in the survey are selected 
from a flux-limited 60 micron a italog (the PSCz; Saun- 
ders et al. 2000), (2) the survey has low and high redshift 
cutoffs at 2; = 0.1 and z = 0.23, and (3) the survey can- 
not include objects close to z = 0.174 where the OH lines 
are redshifted into the strong Galactic HI emission. The 
redshift cutoffs are imposed by radio frequency interfer- 
ence (RFI) encountered at Arecibo above 1510 MHz and 
below 1355 MHz. The typical RMS flux density for a 12 
minute integration is 0.65 mjy, and detections are gen- 
erally made at the 3a level or greater, which is 2 mJy. 
These constraints are well-illustrated by Figure 1, which 
shows redshift cutoffs, the observable candidates, and the 
sensitivity limits in both OH line luminosity and FIR flux. 

Ignoring constraint (1) for now, computation of the OH 
LF can follow the prescription outlined in §3.1 by setting 
the redshift bin Az to span z = 0.1 to z = 0.23 with 
z-min = 0.1. We use na,y = 2.0 mJy and can effectively 
ignore the thin shell of space centered on z — 0.174 since 
its contribution to the total volume is negligible. The sur- 
vey solid angle spans 0° < S < 37° such that fl = 3.78 
steradians (30% of the sky). 

Now fold in the PSCz selection criteria. First, the PSCz 
does not completely cover the 0° < 6 < 37° band. It 
excludes the galactic plane and areas with inadequate or 



12.5 F ^ 




1 ' — ' — • — • — ■ — • — ' — • — • — • — • — ' — • — • — • — 
0.1 0.2 0,3 

z 

Fig. 1. — The Arecibo OH megamaser survey sample, observable 
candidates, and detected OH Megamasers. Shown are the 60 ^im 
luminosities and the redshifts of the complete sample of OHM candi- 
dates extracted from the PSCz (top), the observable subset of these 
candidates which have unambiguous OH line properties {middle), 
and the OH line luminosities of the detected OHMs versus redshift 
(bottom). The dotted vertical lines indicate the 2 = 0.1 cutoff, and 
the solid lines indicate the flux density limits of the PSCz and the 
OHM survey. 



confused lists' coverage (Saunders et al. 2000). The PSCz 
mask excludes 18% of the Arecibo coverage, reducing the 
survey solid angle to f2 = 3.21 steradians (25.5% of the 
sky). The survey volume from 2; = 0.1 to .2 = 0.23 is thus 
0.63 Gpc3. 

Second, the PSCz has a 60 /xm flux density limit of 0.6 
Jy. Hence, the volume available to each object in the sur- 
vey can potentially be limited by this cutoff rather than 
the OH line detectability. The true volume available to a 
given object in the survey is now 

K =min{I/O«,-Kf0Mm}. (11) 

Computation of l/J'O^m fQ^Q-^^g ^i^q continuum prescrip- 
tion outlined in §3.2, and the calculation of the LF uses 
this more restrictive definition of Va- The luminosity bins 
A log L refer to the integrated OH line luminosities of the 
OHMs. Hence, the details of detecting OHMs from both 
surveys are folded into Va, and the luminosity intervals 
in the LF incorporate the information about the OH line 
luminosities. 

Third, the IRAS 60 nm flux measurements require net 
k-corrcctions from to DL^max- The k-corrcctions them- 
selves are derived from spectral energy distribution mod- 
els for star forming galaxies developed by Dale et al. 
(2001). The models depend on the rest- frame FIR color 
log(/6o^m/./ioo^m) which is corrected from the observed 
color to the rest frame color using values tabulated by 
D. Dale (2001, private communication). The typical color 
correction from 2 = 0.15 to 2; = for the OHM host sample 
is 0.05-0.10 (colors get "warmer"). The net k-correction 
for an OHM detected at z = 0.15 which is detectable to 
z = 0.20 with rest frame FIR color of -0.10 is k = 0.98. 
Although the OHM sample spans a wide range of FIR col- 
ors, the net k-corrections vary little from source to source 
in the range z = 0.1-0.23. Objects with no 100 /im detec- 
tion generally have a range of possible net k-corrections 
centered on unity, which is adopted for lack of better in- 
formation. 

Figure 2 shows the OH luminosity, FIR luminosity, and 
redshift distributions of the OHM sample. Also shown are 
the maximum detectable redshift distributions computed 
from the OH detections, the 60 fim detections, and the 
available redshift distribution Za- Note that the Za dis- 
tribution is quite flat but that the maximum detectable 
redshift distribution for OH indicates two populations of 
OHMs: a population which is just detected by the survey 
and indicates a large population of OHMs which would be 
detected by a deeper survey, and a population of "over- 
luminous" OHMs which can be detected in short integra- 
tion times out to z ~ 1, similar to QSOs or radio galaxies. 
Nearly half of the OHM sample falls into this interesting 
latter category. These are the objects which are useful for 
studies of galaxy evolution out to large redshifts. There 
is a hint of this population dichotomy in the plots of Lqh 
versus Lpj^ analyzed in Paper III. 

We combine all of the survey constraints and compute 
an OH LF which is presented in Figure 3. A power-law fit 
to the well-sampled OH LF points gives 

$ ^ {9.8t^','> X 10-«) Loff''*"-'' Mpc--' dex-i (12) 
where Lqh is expressed in solar luminosities. The unifor- 
mity of the sampling in space is checked with the {V/Va) 
test (Schmidt 1968). A uniformly distributed sample be- 
tween and 1 has mean 0.5. Hence, the (V/Va) values 
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Fig. 2. — The OH megamaser sample: Lqh, ^eo^jm, redshift, and maximmn detectable redshift distributions. Panels show; (a) The Lqh 
distribution of the survey OHM detections; (6) The Lsojim distribution; (c) The redshift distribution; (<i) The maximum detectable redshift 
distribution calculated from the OH emission line; (e) The maximum detectable redshift distribution calculated from the 60 /im flux density; 
and (/) The available redshift distribution computed from min{z^^, 2^^^}. These distributions exclude two OHMs at z > 0.23. 



consistent with 0.5 in well-sampled Lqh bins shown in 
Figure 3 indicate a uniformly distributed sample of OH 
mcgamascrs. Both the OH LF and (V/Va) are tabulated 
in Table 3, which includes the number of OHMs available 
in each Lqh bin. 

OH LFs previously computed by Baan (1991) and 
Briggs (1998) show similar properties to the Arecibo OHM 
survey LF, as indicated in Figure 3. Baan's OH LF sam- 
ples a larger range of Lqh, showing a knee at roughly 
10^'^ Lq. The number of OHMs in the Arecibo sample 
contributing to the OH LF below Lqh = W^'^Lq is in- 



adequate to confirm this turnover in the LF. The survey 
sensitivity cutoff at these low line luminosities is severe, 
as seen in Figures 1 and 3. Also noteworthy is the higher 
OHM density found in the Arecibo LF versus Baan's 1991 
result for Lqh > W^Lq. Briggs (1998) derived an OH 
LF from a quadratic OH-FIR relation combined with a 
60/im luminosity function derived by Koranyi & Strauss 
(1997). Briggs computes the OH LF analytically assuming 
an OHM fraction of unity in LIRGs for all L^Qf^m 

(Figure 

3b, dotted line) and numerically for an increasing OHM 
fraction versus ieo/iim (Figure 3b, histogram). Although a 
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Fig. 3. The OH lucgaiiiaser luminosity function. Panels show: (a) The Lqh distribution of the survey OHM detections; (b) The OH 



luminosity function (filled circles) and a power-law fit to the well-sampled data points (<I> oc 



dashed line), the OH LF derived by Baan 



(1991; open circles), and the OH LF computed by Briggs (1998) for a fixed and variable OHM fraction in LIRGs {dotted line and histogram, 
respectively); (c) The redshift distribution; and (d) The average ratio of V to Va which is consistent with 0.5 in well-sampled bins (a uniformly 
distributed sample has (V/Va) = 0.5). These distributions and calculations exclude two OHMs at z > 0.23. 



quadratic OH-FIR relation is not supported by the known 
OHMs (Paper III), Briggs' OH LF follows tlieArecibo OH 
LF data points remarkably closely above Lqh = IO^Xq. 
Briggs obtains a rough power-law slope of —1.15, which is 
inconsistent with the Arecibo result of —0.64 ± 0.2 (Equa- 
tion 12). The inconsistency can be attributed to the steep 
slope at the high luminosity end of the 60/im LF of Ko- 
ranyi & Strauss (1997) compared to the shallower slope 
obtained by Yun, Reddy, & Condon (2001) and Kim & 
Sanders (1998). When translated into an OH LF, the 
steep 60/im slope is significantly lessened by Briggs's use 
of a quadratic OH-FlR relation. The Arecibo power-law 
slope for OHMs is consistent with the power-law LF Kim 
& Sanders (1998) determined for ULIRGs. They found 
$ oc LJ^;^^^°-^ Mpc-3 mag-i over Ljr = IO^^.^qIS^^ 
(the exponent becomes — 0.94± 0.12 when $ is expressed 
in Mpc~^ dex~^). 

A well-determined OH luminosity function forms the 
foundation of any galaxy evolution study which uses OH 
megamasers as luminous radio tracers of mergers, dust- 
enshrouded star formation, or supermassive black hole bi- 



naries. We now use this new OH LF to predict the de- 
tectability and abundance of OHMs available to deep radio 
surveys. 

4. DETECTING OH MEGAMASERS AT HIGH REDSHIFT 

OH megamasers are excellent luminous tracers of merg- 
ing galaxies. They may become a tool for measuring the 
merging history of galaxies across much of cosmic time and 
determine the contribution of merging supermassive black 
holes to the low frequency gravitational wave background 
and to low frequency gravitational wave bursts. They also 
provide an extinction-free tracer of bursts of highly ob- 
scured star formation and may provide an independent 
measure of the star formation history of the universe. Ap- 
plication of OH megamasers to these topics requires (1) 
that they be detectable at moderate to high redshifts, and 
(2) that their sky surface density be high enough for radio 
telescopes to detect at least a few per pointing. 

4.1. Detectability with Current and Planned Facilities 
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Table 3 

The oh LrMixosiTY Ftkctiok 



log Lqjj 


N(OHMs) 


Mpc~^ dex-i 




(V/Va) 


1.6 


1 


(1.5 ± 1.5) X 10- 


7 


0.97± ■•• 


2.0 


2 


(5.4 ±3.9) X 10" 


8 


0.26 ± 0.21 


2.4 


17 


(4.1 ± 1.3) X 10" 


7 


0.62 ±0.16 


2.8 


12 


(1.0 ±0.4) X 10- 


7 


0.48 ±0.16 


3.2 


8 


(8.2 ± 3.3) X 10" 


8 


0.58 ± 0.21 


3.6 


9 


(6.3 ±2.5) X IQ- 


8 


0.49 ±0.19 


4.0 


1 


(3.9 ±3.9) X 10- 


9 


0.83± ••• 



The dctcctability any given OH megamaser at high red- 
shift depends on the strength of the OHM, the sensitivity 
of the instrumentation, and cosmology. Moving beyond 
z '^0.2 requires a careful treatment of cosmology because 
the differences between luminosity distances and volumes 
among the manifold of possible cosmologies become signif- 
icant. We assume for this analysis that Ho = 75 km s~^ 
Mpc-^ Qm = 0.3, and Oa = 0.7. 

Assumptions are also required to translate the observed 
quantity, flux density, to a line luminosity. We assume 
that the integrated flux density can be approximated by 
the product of the peak flux density and an average rest- 
frame width, narrowed by the redshift: 

i'OH — JOH — JOH —7— 7 

1 + z c[l + z) 

where the assumed rest- frame width is Avo = 150 km s~^. 

Figure 4 plots sensitivity thresholds for a 3(7 detection 
of an OHM with luminosity Lqh at redshift z. Included 
are the Arecibo OHM survey detections and sensitivity 



Frequency (MHz) 
400 




AO OHM Survey 
(12 min. 3a) 

Detected OHM 



GMRT (16 hr, 3<j) 

■-■ LOFAR (16 hr, 3a) 
- SKA (12 min, 3a) 



limit, predictions for a 16 hour integration by the Giant 
Metrewavc Radio Telescope (GMRT) at 610 and 327 MHz, 
predictions for a 12 minute integration by the Square Kilo- 
meter Array (SKA) down to 300 MHz, and predictions for 
a 16 hour integration by the Low Frequency Array (LO- 
FAR) below 300 MHz. The GMRT is assumed to span 32 
MHz with 0.12 MHz channels and obtain 0.083 mJy RMS 
noise per beam {Tsys 100 K; Gain = 0.32 K Jy"^).^ The 
SKA sensitivity is computed by scaling up the Arecibo 
OHM survey by a factor of 25 in collecting area but keep- 
ing all other parameters the same. Such a survey would 
have an RMS noise level of 0.027 mJy in 12 minutes of in- 
tegration per 49 kHz channel. The sensitivity of LOFAR 
below 300 MHz is estimated from guidelines provided by 
M. P. van Haarlem'^ to be 0.047 mJy RMS per 0.125 MHz 
channel in 16 hours of integration. 

Note that, RFI environment and available receivers 
aside, Arecibo could detect OH gigamasers out to roughly 
z = 1 in a 12 minute integration. The GMRT shows much 
promise for the detection of Lqh > OHMs at 

z = 1.7 in an integration of 16 hours. In integration times 
of less than an hour, a SKA would be able to detect a large 
fraction of the OHM population out to medium redshifts 
and all OH gigamasers down to its lowest proposed oper- 
ating frequency near 300 MHz. LOFAR would be able to 
detect OH gigamasers in roughly 32 hours of integration 
from z ~ 4.5 back to the reionization epoch if they ex- 
ist. Clearly OH megamasers are detectable at moderate 
redshifts with current facilities and at high redshifts with 
future arrays, but how abundant might they be? 

4.2. The Sky Density of OH Megamasers 

The OH megamaser luminosity function (OH LF) de- 
rived in §3 can predict the sky density of detectable OHMs 
as a function of instrument sensitivity, bandpass, and red- 
shift. A useful function in terms of observational param- 
eters would be the number of OHMs detected per square 
degree on the sky per MHz bandpass searched. This can 
be expressed as an integral of the OH LF over the range 
of detectable Lqh- 

flog LoH,max 

d log Lqh 



dN 
dfi df 



J 

Jlot 

L 



/iogiofr.™ir.(z) dVLdvd log LoH 
dN dV 



d log LoH 



Fig. 4. — Detectability of OH megamasers. Contours show the 
sensitivity required to detect an OHM with luminosity Lqh at red- 
shift z. Included are the results and sensitivity of the Arecibo OHM 
survey and predictions for the GMRT, the SKA, and LOFAR. 



/ dV d log Lqh dQdu 

^ See http:/ /www.ncra.tifr.rcs.in/ncra_hpagc/gmrt/gmrt_spec.html. 
^ See http;//www. lofar.org/science/index. html. 
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dV 



^{Loh) d log LoH- 



(14) 



Frequency {MHz) 
600 400 



dil dv 

Recall that the OH LF is the number of OHMs with lu- 
minosity LoH per Mpc^ per logarithmic interval in Loh, 
expressed as 

$(W) = 6io//- (15) 
Hence, the integral in Equation 14 reduces to a simple 
form 

^{Loh) dlogLoH = [LoH,max ' LoH,mM] ■ 

The volume clement per unit solid angle per unit frequency 
can be translated into the usual panoply of cosmological 
parameters (Weinberg 1972): 
dV _ dz dV 
dQdu dv d^dz 

_ dV {1 + zf 
dfl dz Vo 

= , . (17) 

Hence, the final result for the sky density of OHMs in a 
flat spacetime has the form 

dN cD'-, b 



dO. dv HoVo ^Jil + zyVlM + ^A « In 10 



^ \_LoH,max ^ H ,min{^)\ ■ (18) 

The minimum detectable Lqh depends on the sensitivity 
of observations and the redshift, as shown in Figure 4. The 
maximum Lqh is assumed to be 1O^ '*L0 which is a fac- 
tor of two larger than the luminosities of the known OH 
gigamascrs. The ofFccts of relaxing this conservative as- 
sumption are discussed below. When LoH,min = LoH,max 
there are no OHMs left to detect. 

Finally, we fold in an evolution parameter m which 
scales the sky density of OHMs as (1 + z)"\ Assuming 
Ho = 75 km's^i Mpc-\ I^m = 0.3, and f^A = 0.7, and 
folding in the OH LF parameters a and b from §3, we 
obtain the sky density of OHMs for any sensitivity level 
as a fimction of redshift. The detectable sky density of 
OHMs versus redshift is plotted in Figure 5 for a 3a OH 
line detection at 0.25 mJy. Several evolution scenarios 
are plotted: m = (no evolution), m = 4, and m = 8. 
We assume an arbitrary turnover in the evolution factor 
at z — 2.2, after which the number density of OHMs is 
constant. This turnover in the merger rate corresponds 
roughly to the turnover in the QSO luminosity function. 
The point with error bars is a prediction of the sky density 
of OHMs derived from the space density of submillime- 
tcr detections (which arc ostensibly ULIRGs) computed 
by Barger, Cowie, & Richards (2000) assuming that half 
of the submillimeter galaxies produce detectable OHMs. 
Bands currently accessible to the GMRT (610 and 327 
MHz) in roughly 16 hours of integration are indicated in 
bold. The proposed Square Kilometer Array might reach 
the same sensitivity level in a few minutes of integration. 
If the submillimeter predictions are to be believed, very 
strong evolution is favored and the GMRT could detect 
dozens of OHMs per square degree in its 32 MHz band- 
pass in a reasonable integration time. This is a very ex- 
citing prospect and can indicate which models of galaxy 
evolution are favored, as discussed below. 



2 - 








Fig. 5. — Sky density of detectable OH megamasers at 0.25 mJy. 
The expected detections for several galaxy merger histories are la- 
beled by the rate of increase in the merger rate as (1 + z)"^. The 
point with error bars is the prediction for OHM detections based on 
the space density of submillimeter detections computed by Barger, 
Cowie, & Richards (2000). The turnover in merger rate from in- 
creasing to constant at z = 2.2 is indicated by the dotted vertical 
line. The GMRT (currently accessible bands are indicated in bold) 
can reach the RMS noise of 0.083 mJy in each 0.12 MHz channel 
in roughly 16 hours of integration. The proposed SKA might reach 
the same noise level in a few minutes of integration while LOFAR 
would require several hours. 



Briggs (1998) makes predictions of the sky density of 
OHMs using the FIR LF of LIRGs of Koranyi k Strauss 
(1997), a quadratic OH-FIR relation, and an optimistic 
OHM fraction as a function of igopm derived from Baan 
(1989) and Baan et al. (1992). Briggs predicts roughly 16 
OHMs at 2; = 2 per square degree per 50 MHz bandpass for 
m = 4.5 and a sensitivity level of 0.2 mJy. The predictions 
of OHM detections based on the Arecibo survey OH LF are 
not as rosy; only about 6 OHMs would be detected for m = 
4.5 at 0.2 mJy &t z = 2. The main discrepancy between 
the two predictions lies in Briggs' use of a quadratic OH- 
FIR relation and a high OHM fraction in LIRGs compared 
to the fraction detected in the Arecibo OHM survey. 

If the assumptions built into these detection predictions 
are changed, the results change only slightly below z ^ 
2. For example, if the OH LF is extended to Lqh = 
IO^Lq, there is no significant effect on the detection rates 
below z ~ 3 because the most luminous OHMs are rare. 
The detection rates at the highest redshifts will change 
significantly because only the tail of the OH LF is detected, 
but predictions beyond z = 2 are wild speculations with 
no supporting data. The choice of cosmology also docs 
not alter the detection rates dramatically. Changing from 
an Qa = to an Oa = 0.7 universe amounts to roughly 
a factor of 1.5 increase in OHM detections at z < 2. The 
effect again becomes more significant at higher redshifts. 
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4.3. The Merging History of Galaxies 

A property as basic as the number density of galaxies 
as a function of rcdshift is not well known. The number 
density of galaxies depends of course on the merging his- 
tory of galaxies, which is an essential ingredient in theories 
of structure formation and galaxy evolution. There have 
been a number of studies which attempt to measure the 
merger rate or merger fraction of galaxies versus redshift, 
and most of them parameterize the evolution of mergers 
with a factor (1 + z)"^. A Hubble Space Telescope survey 
identifying optical mergers and close pairs of galaxies in- 
dicates m = 3.4 ± 0.6 from 2; = to ^ = 1 (Le Fevre 
et al. 2000), whereas the 1 Jy ULIRG survey indicates an 
exponent of m = 7.6 ± 3.2 (Kim & Sanders 1998) which 
is similar to the evolution of bright QSOs with redshift as 
m - 6 up to z - 2.5 (Briggs 1998; Hewett et al. 1993; 
Schmidt et al. 1995). Semi-analytic models of hierarchical 
galaxy formation by Kauffmann & Haehnelt (2000) indi- 
cate an evolution of the number density of gas-rich merg- 
ers from 2; = to 2; = 2 roughly following m = 4. The 
sampling of work on this topic listed here is by no means 
complete. In general, most stiidies indicate a strong in- 
crease of the merging rate with redshift, but it is likely 
to be a strong function of the total dark halo mass (for 
which bolometric luminosity is a proxy) and the mass ra- 
tio of the merging pair (Khochfar & Burkert 2001). A 
luminosity dependence is seen in the evolution of quasars 
and LIRGs, with the most luminous objects showing the 
strongest evolution (Schmidt et al. 1995; Kim & Sanders 
1998). 

Although a deep survey for OH megamasers may not 
be translatable into a merger rate, it does discriminate 
between various galaxy evolution scenarios. For example, 
there are two orders of magnitude difference between the 
OHM detections expected from m = 8 versus m = 4 at 
z ^ 2. Even a lack of detections in an adequately deep 
OHM survey would provide new constraints on the evolu- 
tion of mergers with redshift and would suggest a revision 
of the conventional interpretation of submm galaxies as 
high-redshift ULIRGs. If, on the other hand, OHMs were 
detected in a deep survey, this would not only indicate the 
evolution of merging but also provide an extinction-free 
redshift determination for submillimeter galaxies which 
have proven to be extremely difficult to identify and ob- 
serve at optical wavelengths or even IR wavelengths (e.g. — 
Small et al. 2000; Ivison et al. 2000). 

4.4. The Star Formation History of the Universe 

Merging systems traced by OH megamasers are in the 
throes of extreme starbursts, but much of the star forma- 
tion is highly obscured by dust. The history of obscured 
star formation is difficult to determine with current in- 
struments which detect the reprocessed light directly at 
FIR or submm wavelengths due to poor sensitivity or the 
difficulty of determining redshifts, respectively. OH mega- 
masers are a promising solution to these limitations be- 
cause they are luminous, are unattenuated by dust, and 
provide redshifts. As Townsend et al. (2001) point out, 
OH megamasers may be used to determine the nature and 
evolution of submm galaxies and to indicate their rele- 
vance to the star formation history of the universe. If the 
connections between OHMs and their hosts, particularly 



between Lqh and LpiB, (or at least the OH LF and the 
FIR LF), are valid at moderate to high redshifts, then sur- 
veys for OHMs can be used as an independent measure of 
the star formation history across a large span of the age 
of the universe. 

4.5. Constraining The Gravitational Wave Background 

One final potential application of OHMs addresses the 
end result of massive mergers: the formation and coales- 
cence of binary supermassive black holes. Supposing that 
each galaxy in a merger contains a supermassive black 
hole, the rapid coalescence of nuclei due to dynamical fric- 
tion will produce a binary supermassive black hole which 
will continue to decay until a final coalescence event which 
will produce (among other things) a burst of gravitational 
waves. Bursts from supermassive black holes are likely to 
be the major source of 10~^-10° Hz gravitational waves 
(Haehnelt 1994) which may someday be detectable by 
the Laser Interferometer Space Antenna (LISA) or long- 
duration pulsar timing. The merging rate of galaxies, as 
well as the masses of the black holes involved, determine 
the event rate detectable by LISA. The integrated merg- 
ing history of galaxies determine the noise levels produced 
in pulsar timing, and would provide a "foreground" to 
any cosmological background of gravitational waves pro- 
duced during the inflationary epoch (D. Backer 2001, pri- 
vate communication). Clearly, getting some handle of the 
event rate of supermassive black hole mergers would pro- 
vide much needed constraints and thresholds for the diffi- 
cult work of detecting gravitational waves. 

5. SUMMARY 

The OH luminosity function constructed from a sample 
of 50 OH megamasers detectcsd by the Arecibo OH mega- 
maser survey indicates a power-law falloff with increasing 
OH luminosity 

$ = [Q.Stf.^^ X 10-6) L50^64±o.2i 

valid for 2.2 < log Lqh < 3.8 (expressed in Lq) and 
0.1 < z < 0.23. The OH LF is used to predict the areal 
density of detectable OHMs at arbitrary redshift for a 
manifold of galaxy merger evolution scenarios parameter- 
ized by (1 -I- 2)™ where < m < 8. For reasonable choices 
of m, an "OH Deep Field" obtained with the Giant Aletre- 
wave Radio Telescope at 610 MHz {z = 1.73) may detect 
dozens of OHMs per square degree in a reasonable integra- 
tion time. A lack of detections in a sufficiently deep field 
would also significantly constrain the evolution of merging 
and exclude the most extreme evolution scenarios. 
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